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[ij  The  Navy  Coastal  Ocean  Mode!  (NCOM)  is  a  frec-surface,  primitive-equation 
model  that  is  under  development  at  the  Naval  Research  Laboratory  (NRL).  The 
NCOM-bascd  model  of  the  Monterey  Bay  area  is  evaluated  during  a  series  of  upwelling 
and  relaxation  wind  events  in  August- September  of  2000.  The  model  receives  open 
boundary  conditions  from  a  regional  NCOM  implementation  of  the  California  Current 
System  and  surface  fluxes  from  the  Navy  Coupled  Ocean/Atmosphcrc  Mcsoscale 
Prediction  System  (COAMPS™)(COAMPS  is  a  registered  trademark  of  the  Naval 
Research  Laboratory).  Issues  investigated  in  this  study  are:  NCOM-bascd  model 
simulations  of  upwelling  and  relaxation  events,  coupling  to  CO  AM  PS,  use  of  sigma 
versus  hybrid  (sigma-z)  vertical  grids,  and  coupling  with  a  largcr-scalc  model  on  the  open 
boundaries.  The  NCOM  simulations  were  able  to  reproduce  the  observed  sequence  of 
the  upwelling  and  relaxation  events,  which  can  be  attributed,  in  part,  to  the  good 
agreement  between  the  observed  and  COAMPS  winds.  Comparisons  with  the  mooring 
observations  show  that  COAMPS  overestimates  shortwave  radiation  values,  which  makes 
the  NCOM  modeled  SSTs  too  warm  in  comparison  with  observations.  The  NCOM 
runs  forced  with  different  resolution  atmospheric  forcing  (3  versus  9  km)  do  not  show 
significant  differences  in  the  predicted  SSTs  and  mixed-layer  depths  at  the  mooring 
locations.  At  the  same  time,  during  the  extended  upwelling  event,  the  model  runs  forced 
with  3  and  9  km  resolution  COAMPS  fields  show  differences  in  the  surface  circulation 
patterns,  which  arc  the  most  distinct  in  the  southern  portion  of  the  model  domain. 

The  model  run  with  9-km  forcing  develops  a  northward  flow  along  the  coast,  which  is  not 
present  in  the  run  with  3 -km  forcing  and  in  observations  (for  example,  HF 
radar-derived  radials).  Comparison  of  the  wind  patterns  of  the  3-  and  9-km  products  shows 
a  weakening  of  the  9-km  wind  stress  along  the  southern  coast  of  the  NCOM  model 
domain,  which  is  responsible  for  the  development  of  the  artificial  northward  flow  in  the 
NCOM  run  with  9-km  forcing. 
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1*  Introduction 

[2]  A  multi-institution,  multidisciplinary  field  experiment 
was  conducted  in  the  Monterey  Bay  area  during  August- 
September  of  2000  [Ramp  et  al .,  2005].  The  experiment 
represented  a  joint  effort  of  the  Autonomous  Ocean  Sam¬ 
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pling  Network  (AOSN  I),  the  Monterey  Bay  Aquarium 
Research  Institute  (MBARI)  Upper- Water  column  Science 
Experiment  (MUSE),  and  the  National  Oceanic  Partnership 
Program  (NQPP)  Innovative  Coastal-Ocean  Observing  Net¬ 
work  (ICON)  projects.  The  experiment  objectives  ranged 
from  a  study  and  modeling  of  biophysical  processes  of  the 
upwelling  centers  to  the  development  of  components  of  a 
real-time  observational  and  modeling  system  for  the  Central 
California, 

[3]  The  observational  network  consisted  of  a  combination 
of  airborne,  in  situ,  AUVs,  and  shore-based  HF  radar 
observations,  and  the  observational  period  (15  August  - 
2  September)  included  an  extended  upwelling  event  during 
15-21  August  and  an  extended  wind  relaxation  event  from 
28  August  through  2  September.  The  numerical  ocean 
model  was  based  on  the  Princeton  Ocean  Model  [POM, 
Blumherg  and  Mellor,  1987]  of  the  Central  California  Coast 
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developed  under  the  ICON  project  [Paduan  and  Shulman, 
2004;  Shulman  et  at.  ,  2002]. 

[4]  As  noted  in  Ramp  et  aL  [2005],  the  POM  ICON 
model  results  were  very  consistent  with  observations  during 
the  upwelhng  event  but  differed  significantly  during  the 
extended  relaxation  period.  During  the  relaxation  period, 
strong  onshore  translation  of  warmer  California  Current 
water  combined  with  local  surface  heating  inside  the  Bay 
resulted  in  extraordinary  surface  warming  within  the  bay 
during  30-31  August.  It  was  hypothesized  by  Ramp  et  al 
[2005]  that  the  model’s  difficulty  in  reproducing  the  extreme 
heating  event  was  likely  due  to  inaccuracies  in  both  the 
wind  stress  and  surface  heat  fluxes  from  the  9-km 
COAMPS  and  insufficient  vertical  resolution  of  the  upper 
ocean  in  the  model  itself. 

[5]  The  POM  ICON  model  [used  by  Ramp  et  aL ,  2005]  is 
a  three-dimensional,  free-surface  model  with  I  -4  km  hori¬ 
zontal  resolution  and  thirty-vertical  sigma  levels  [Shulman 
et  al .,  2002].  The  model  was  forced  with  9-km  resolution 
wind  stresses  and  heat  fluxes  from  the  Navy  Coupled  Ocean 
and  Atmospheric  Mesoscale  Prediction  System  (COAMPS™) 
reanalysis  for  the  US  west  coast  [Hoduret  al ,  2002;  Kindle  et 
al,  2002)  and  by  the  regional  Pacific  West  Coast  (PWC) 
model  [Haidvogel  etal. ,  2000;  Rock  ford  and  Shulman,  2000] 
on  the  seaward  boundaries.  The  model  assimilated  surface 
currents  from  HF  radar  (CODAR)  sites  around  the  Monterey 
Bay  area  [Paduan  and  Shulman,  2004]. 

[t>]  In  the  present  paper,  the  Navy  Coastal  Ocean  Model 
(NCOM)  [Rhodes  et  al.,  2002;  Martin,  2000]  together  with 
two  new  high-resolution  COAMPS  surface  flux  fields  are 
utilized  in  order  to  evaluate  the  NCOM  model's  predictive 
capabilities  and  to  examine  the  sensitivity  of  the  simulated 
circulation  patterns  to  the  model  formulation  and  surface 
forcing  functions.  The  NCOM  model  is  based  on  POM,  but 
has  options  for  using  different  turbulence  closure,  advec- 
tion,  and  open  boundary  condition  schemes,  etc.  The 
NCOM  in  Monterey  Bay  (noted  belowf  as  NCOM  ICON 
model)  is  set  up  on  the  same  grid  as  the  POM  ICON.  Since 
August  of  2000,  a  number  of  modifications  were  made  to 
the  CO  AM  PS,  most  important  differences  are  the  addition 
of  a  high-resolution  3-km  nest  and  improvements  to  the 
surface  flux  scheme  and  microphysics  parameterization. 

[7]  Models  can  provide  us  with  additional  information 
about  circulation  patterns  and  water  mass  properties  and 
formations  which  can  not  be  revealed  by  limited  observa¬ 
tions  alone.  The  circulation  patterns  of  the  POM  ICON 
model  simulations  [Ramp  et  al,  2005]  revealed  distinct 
differences  in  the  circulation  and  water  mass  properties  of 
Monterey  Bay  during  the  upwrelling  and  relaxation  events. 
In  the  present  paper,  our  objective  is  to  continue  the  analysis 
of  the  circulation  patterns  and  water  mass  properties  pre¬ 
dicted  by  the  models  of  Monterey  Bay. 

[s]  The  paper  has  the  following  structure:  In  section  2, 
we  describe  the  hierarchy  of  different  resolution  NCOM- 
based  models  in  the  Pacific  Ocean  used  in  this  study.  A 
description  of  the  atmospheric  forcing  is  presented  m 
section  3,  Section  4  provides  a  brief  description  of  the 
physical  conditions  in  the  Bay  during  2  weeks  of  August 
2000  as  well  as  an  evaluation  of  COAMPS  predicted 
surface  fluxes.  Section  5  presents  an  evaluation  and  analysis 
of  the  model  results.  This  section  is  split  into  four  sub¬ 
sections,  which  focus  on  the  model  sensitivity  to  the 


resolution  of  the  atmospheric  forcing  and  to  the  vertical 
coordinate  system  and  simulations  during  the  extended 
relaxation  and  upwelling  events.  Section  6  is  devoted  to 
discussion  and  conclusions. 

2.  Description  of  Hierarchy  of  NCOM-Based 
Different  Resolution  Models  in  the  Pacific  Ocean 

[9}  NCOM  is  a  primitive -equal ion,  three-dimensional, 
hydrostatic  model  [Rhodes  et  al. ,  2002;  Martin ,  2000]. 
The  model  is  based  on  POM,  but  has  options  for  using 
different  turbulence  closure,  advection,  and  open  boundary 
condition  schemes,  etc.  The  model  includes  routines  for 
direct  reading  and  processing  of  atmospheric  COAMPS 
fields.  NCOM  allows  for  the  use  of  a  sigma  or  hybrid 
sigma/z-level  vertical  coordinate  systems.  The  code  is  par¬ 
allelized  and  runs  on  many  different  computer  platforms. 

[  jo]  To  accommodate  the  wide  range  of  horizontal  scales 
of  physical  and  biological  oceanic  processes,  the  Naval 
Research  Laboratory  (NRL)  uses  a  hierarchy  of  different 
resolution,  NCOM -based  ocean  models.  The  NCOM  global 
model  [Rhodes  et  a l 2002]  is  run  in  real-time  and  uses  a 
hybrid  vertical  coordinate  system  with  19  sigma  layers  from 
the  surface  down  to  138-m  depth  and  21  z-levels  between 
138  m  and  the  bottom.  Global  NCOM  has  1/8°  horizontal 
resolution,  assimilates  three-dimensional  temperature  and 
salinity  observations  derived  from  the  Modular  Ocean  Data 
Assimilation  System  [MOD AS;  Fox  et .  al .,  2002],  and  uses 
atmospheric  forcing  from  the  Navy  Global  Atmospheric 
Prediction  System  [NOGAPS;  Rosmond  et  al. ,  2002]. 

[u]  Global  NCOM  provides  boundary  conditions  for  an 
NCOM-based  regional  model  of  the  California  Current 
[NCOM  CCS,  Shulman  et  al .,  2004].  The  NCOM  CCS  has 
a  horizontal  resolution  of  about  9  km  and  is  used  with  two 
different  vertical  coordinate  systems:  (l)  a  purely  sigma- 
coordinate  grid  with  30  layers,  and  (2)  the  same  40-layer, 
hybrid  grid  used  in  the  global  model.  In  the  present  study, 
the  NCOM  CCS  with  the  hybrid  (40-layer)  vertical  coordi¬ 
nate  system  is  used.  The  NCOM  CCS  model  is  forced  with 
atmospheric  products  derived  from  the  COAMPS  predic¬ 
tions  [Kindle  et  aL,  2002].  As  in  the  global  model,  the 
NCOM  CCS  assimilates  three-dimensional  temperature  and 
salinity  observations  derived  from  the  MGDAS. 

[12]  NCOM  CCS  provides  boundary  values  to  the  high- 
resolution  NCOM  ICON  model  of  the  Monterey  Bay  Area, 
The  NCOM  ICON  runs  on  an  orthogonal,  curvilinear  grid, 
which  has  a  variable  horizontal  resolution  ranging  from  !  to 
4  km,  with  finer  resolution  in  and  around  Monterey  Bay 
(Figure  1).  The  NCOM  ICON  is  set  up  on  the  same 
horizontal  grid  as  the  POM  ICON  [Shulman  et  al .,  2002], 
The  model  uses  the  Mellor-Yamada  level  2.5  turbulence 
closure  scheme  [ Martin ,  2000].  In  this  study,  the  NCOM 
ICON  is  forced  with  atmospheric  fluxes  from  COAMPS  (no 
bulk  formulations  for  heat  fluxes  or  relaxation  terms  are 
used,  see  section  3).  The  MODAS  derived  temperature  and 
salinity  are  not  assimilated  into  the  fine-resolution  NCOM 
ICON  model, 

3.  Atmospheric  Forcing 

[13]  The  atmospheric  fields  from  two  COAMPS  versions 
(COAMPS  v2  and  CO  AM  PS  v3)  are  used  in  this  study.  A 
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Figure  1.  The  Monterey  Bay  model  domain  with  locations 
of  moorings  Ml  and  M2. 

number  of  modifications  were  made  to  the  physical  param- 
eterizations  during  a  5-year  period  that  encompassed  the 
advancement  ofCOAMP$v2  to  COAMPSv3.  The  two  most 
important  differences  between  the  two  versions  of  the 
COAMPS  are  related  to  the  surface  flux  and  microphysics 
parameterization  In  particular,  the  modified  surface  flux 
scheme  improves  the  wind  stress  calculation  as  compared 
with  observed  values  [Wang  et  aLy  2002].  Consequently,  the 
new  scheme  significantly  improves  the  near- surface  wind 
prediction  as  the  wind-speed  bias  is  reduced  by  half  during 
the  daytime  heating  cycle  over  land  areas.  Additionally, 
improvements  to  the  surface  energy  budget  over  land 
(including  the  soil  moisture  representation  and  deep  soil 
parameterization)  were  introduced.  The  microphysics  pa¬ 
rameterization  is  formulated  following  Rutledge  and  Hobbs 
[1983].  A  number  of  improvements  were  made  to  the 
microphysics  including  a  new  representation  of  the  ice 
phase  microphysics  and  a  reformulation  of  the  entire 
scheme  in  order  to  achieve  total  water  conservation. 

[J4]  COAMPS  assimilates  atmospheric  observations  from 
radiosondes,  aircraft,  satellite,  and  ships  [Hodur  et  al , 
2002],  The  COAMPS  SST  analysis  is  performed  directly 
on  the  particular  nest  grid  and  includes  assimilation  of 
observations  from  ships,  buoys,  satellites  (for  example, 
Multi-Channel  Sea  Surface  Temperature  (MCSST)).  Atmo¬ 
spheric,  as  well  as  oceanic  observations  from  moorings  M 1 
and  M2  (Figure  1)  were  not  assimilated  into  the  COAMPS 
simulations, 

[15]  The  3- and  9-km  resolution  nests  from  the  COAMPSv3 
version  will  be  referred  to  as  "3km,  Cv3"  and  "9km,  Cv3*\ 
respectively.  The  9-km  nest  from  the  COAMPS v2  will  be 
referred  to  as  "9km,  Cv2”.  All  COAMPS  products  provide 
atmospheric  fluxes  with  hourly  frequency  and  the  NCOM 
interpolates  them  to  the  current  model  lime. 

4,  Brief  Description  of  Observed  Physical 
Conditions  and  COAMPS  Predictions  During 
15  August -2  September  2000 

[ig]  A  plot  of  observed  wind  velocities  at  moorings  Ml 
and  M2  (Figure  1)  is  presented  in  Figure  2a.  The  plot 
indicates  that  during  15=21  August  there  was  a  strong 
upwelling  favorable  event,  which  was  followed  by  a  brief 


relaxation  event  of  22=23  August.  Upwelling  favorable 
w-inds  resumed  during  24=27  August.  This  brief  upwelling 
event  was  followed  by  an  extended  relaxation  event  of 
28  August- 2  September  2000  [see  also,  Ramp  et  a/.,  2005 
and  Fitzwater  et  al .,  2003].  Our  major  focus  will  be  on  the 
model's  ability  to  reproduce  the  extended  relaxation  event 
of  28  August-2  September  and  the  upwelling  event  of 
17=21  August 

[17]  During  upwelling,  cold  offshore  flowing  filaments 
can  be  seen  at  two  upwelling  centers:  north  of  the  Monterey 
Bay  around  Pt.  Ano  Nuevo,  and  south  at  Pi,  Sur  (Figure  1). 
In  between  these  upwelling  centers,  there  often  exists  a 
warm  mesoscale  eddy.  The  position  of  the  eddy  is  observed 
to  migrate  onshore  and  offshore.  The  upwelling  event  of 
17-21  August  was  characterized  by  the  appearance  of  cold, 
salty  water  at  Point  Ano  Nuevo  (Figure  1)  which  subse¬ 
quently  spread  southward  across  the  mouth  of  the  bay  as  the 
winds  continued.  The  circulation  in  the  bay  itself  was 
cyclonic  with  very  warm  water  observed  in  the  northeast 
comer  of  the  bay.  During  the  wind  relaxation  event,  the 
surface  current  and  temperature  response  was  dominated  by 
the  onshore  translation  of  the  offshore  eddy  and  by  local 
surface  heating  in  the  bay  itself.  Surface  currents  in  the  bay 
remained  cyclonic. 

[j&]  Plots  of  wind  velocities  from  the  different  COAMPS 
products  at  stations  M 1  and  M2  are  presented  in  Figures  2b, 
2c,  and  2d.  Figures  2b,  2c,  and  2d  also  show  the  complex 
correlations  (a)  and  angular  displacemenLs  (theta)  [estimated 
according  to  Kundu,  1976]  between  the  observed  wind 
velocity  and  the  COAMPS -predicted  wind  velocities  for 
the  three  products  considered  here.  The  angular  displace¬ 
ment  gives  the  average  counterclockwise  angle  of  the 
COAMPS  wind  velocity  with  respect  to  the  observed  wind 
velocity.  All  the  CO  AM  PS  products  show  a  correlation 
higher  than  0.68  with  the  observed  wind  velocity  at  Ml  and 
M2  and  small  angles  of  displacement.  The  correlation  0.68 
is  significant  (different  from  zero)  if  the  number  of  degrees 
of  freedom  is  larger  than  or  equal  to  5  for  a  significance 
level  of  0.05  (see  for  example.  Table  XI  of  Hogg  and  Tanis 
[1993]).  The  correlations  on  Figure  2  were  calculated 
betw  een  the  hourly  COAMPS  and  observed  wind  velocities. 
Therefore  around  432  samples  were  used  to  compute  the 
correlations.  It  is  clear  that  the  number  of  degrees  of 
freedom  is  smaller  than  432,  but  we  can  assume  that  it  is 
larger  than  5  for  the  following  reasons:  two  upwelling  and 
two  relaxation  events  are  observed  over  the  considered  time 
frame;  there  is  a  strong  diurnal  variability  in  the  atmospheric 
conditions  in  the  area,  and  correlations  over  an  18-day  time 
frame  were  estimated.  At  the  same  time,  we  cannot  con¬ 
clude  that  correlations  on  Figure  2  are  significantly  different 
due  to  uncertainty  about  the  actual  number  of  degrees  of 
freedom  and  small  differences  between  the  correlations  for 
the  three  considered  CO  AM  PS  products. 

[19]  There  is  remarkable  agreement  between  CO  AMPS 
predicted  winds  and  observations  in  the  sequence  and 
extent  of  each  upwelling -relaxation  event  observed  during 
15  August-2  September.  However,  at  Ml  there  is  an 
overestimation  of  wind  velocity  during  the  first  upwelling 
for  the  “3km,Cv3^  model  run.  The  nearshore  winds  in 
"3km, Cv3"  are  probably  amplified  in  comparison  to 
"9km,Cv3"  because  of  the  better  resolution  of  the  terrain 
and  coastline  features.  The  coastal  topography  is  steeper  in 
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August  -  September  2000  HOURLY  WINDS 
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9km,  Cv2  at  M2 
a=0.83,  theta=8,54 


9km, Cv2  at  Ml 
a=0.73,  theta=5.54 


Figure  2*  Observed  and  COAMPS  model  predicted  wind  velocities  at  Ml  and  M2:  (a)  observations: 
(b)  COAMPS  version  3  with  3-km  resolution  (3km,Cv3);  (c)  COAMPS  version  3-  with  9-km  resolution 
(9km, Cv3);  (d)  COAMPS  version  2  with  9-km  resolution  (9km, Cv2);  Notations  “Up,"  and  “Re."  mean 
upwelling  and  relaxation. 


the  3-km  mesh  and  can  produce  wakes  and  wakes  flow 
regions  with  sharper  gradients.  This  is  particularly  the  case 
in  the  strong  northwesterly  flow  regime  in  which  the  coastal 
terrain,  such  as  the  Santa  Cruz  Mountains,  can  force  a  lee- 
side  w  ake.  The  acceleration  as  a  result  of  the  boundary  layer 
flow  around  the  headlands  can  be  larger  as  well  due  to  the 
steeper  and  better  resolved  terrain. 

[20]  Plots  of  daily  averaged,  photosynthetically  available 
radiation  (PAR)  observed  at  MI  and  M2  and  estimated  from 
the  COAMPS  short  wave  radiation  (SWR)  fluxes  are  pre¬ 
sented  in  Figure  3a  for  Ml  and  Figure  4a  for  M2.  The 
observed  PAR  was  measured  by  the  Biospherical  PRR-620 
spec  trorad  tome  ter  mounted  on  moorings  approximately  3  m 


above  the  water  surface  [Chavez  el  ai ,  2000],  PAR  from 
COAMPS  predictions  was  estimated  as  45%  of  the 
CO  AM  PS-predicted  SWR  flux  [Sirutfon  and  Chavez , 
2004].  Anomalies  of  the  observed  and  model -predicted 
PAR  values  (deviations  from  the  mean  value  over  the 
considered  time  period)  are  plotted.  Observed  and  model 
predicted  mean  values  of  PAR  are  shown  in  the  upper  left 
comer  of  each  subfigure  of  Figures  3a  and  4a. 

[21]  From  the  observed  and  model -predicted  mean  values 
presented  in  Figures  3a  and  4a,  it  is  clear  that  there  is  an 
overestimation  of  the  SWR  in  the  CO  AMPS  predictions  for 
all  the  CO  AM  PS  products  considered  here,  especially  for 
the  “3km, Cv3"  and  “9km, Cv3"  products.  The  excessive 
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Figure  3,  (a)  Observed  and  the  CO  AMPS -derived  PAR  anomalies  at  M 1  location  from  three  CO  AM  PS 
products:  3km,Cv3;  9km,Cv3;  and  9km,Cv2  (see  captions  of  the  Figure  2  for  explanations  of  notations), 
(b)  Observed  and  the  model  predicted  SSTs  from  the  NCOM  ICON  runs  1-3  (see  Table  1)  forced  with 
corresponding  CO  AM  PS  products  shown  on  the  left  panels. 


Anomaly  of  Daily  Mean  Anomaly  of  Daily  Mean 

Surface  PAR  at  M2  (W/m2)  Surface  Temperature  at  M2  (°C) 
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Figure  4.  Same  as  Figure  3  but  at  mooring  M2. 
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Tablet.  Description  of  NCOM  Runs 


Runs 

COAM  PS  Version 

COAMPS  Resolution 

SWR  Correction 

Vertical  Grid 

Run  1 

COAMPSv3 

3  km 

None 

30  sigma  layers 

R  un  2 

COAMPSv3 

9  km 

None 

30  sigma  layers 

Run  3 

COAMPSv2 

9  km 

None 

30  sigma  layers 

Run  4 

COAMPSv3 

3  km 

Reduced  by  75  W/m2 

30  sigma  layers 

Run  5 

COAMPSv3 

3  km 

Reduced  by  a  factor  of  1 .4 

30  sigma  layers 

Run  6 

COAMPSv3 

9  km 

Reduced  by  a  factor  of  1 ,4 

30  sigma  layers 

Run  7 

COAMPSv3 

3  km 

Reduced  by  a  factor  of  1 .4 

Hybrid,  40  layers  (19  sigma,  2l-z) 

Run  & 

COAMPSv3 

9  km 

Reduced  by  a  factor  of  1 .4 

Hybrid,  40  layers  (19  sigma,  2  \~z) 

SWR  is  likely  related  to  a  modeling  of  low-level  clouds, 
which  are  very  extensive  during  summer  time  in  the 
Monterey  Bay  area.  Accurate  prediction  of  low-level  clouds 
is  a  very  challenging  problem.  We  believe  that  the 
COAMPS  model  probably  underestimates  low-level  clouds 
in  the  Monterey  Bay  area. 

[22]  Analysis  of  the  anomalies  presented  in  Figures  3  a 
and  4a  exhibits  a  good  correspondence  between  the 
COAM PS-predicted  and  observed  PAR,  especially  for  the 
“3km, Cv3"  and  “91tm*Cv3”  products.  The  older  version  of 
the  COAM  PS  “9  km,Cv2”  run  completely  missed  the  local 
increase  in  SWR  observed  during  the  relaxation  event  of 
27-31  August  (see  Figure  3a,  bottom  panel).  As  shown 
below,  this  missed  event  in  the  predicted  SWR  was  one  of 
the  reasons  why  the  POM  ICON  model  was  not  able  to 
reproduce  the  extraordinary  surface  heating  within  the  bay 
during  30-31  August  [Ramp  et  al.t  2005). 

5,  Results 

[23]  Table  l  describes  attributes  of  the  NCOM  ICON  runs 
evaluated  in  the  paper  All  the  considered  model  runs  were 
initialized  on  August  2  of  2000  by  using  fields  ftom  the 
same  NCOM  CCS  model,  which  was  forced  with  the 
“9km,Cv3  ^atmospheric  forcing.  NCOM  ICON  results  during 
the  1 5  August -2  September  are  evaluated, 

5,1,  Sensitivity  to  the  Resolution  of  the 
Atmospheric  Forcing 

[24]  First,  we  compare  the  following  NCOM  ICON  runs 
(see  Table  1):  run  1,  the  model  is  forced  with  “3km,Cv3"'; 
run  2,  the  model  is  forced  with  “9km,Cv3,,;  and  run  3,  the 
model  is  forced  with  “9km,Cv2". 

[25]  Figures  3b  and  4b  present  observed  and  model- 
predicted  SST  anomalies  (deviations  from  the  corresponding 
means  over  considered  time  frame)  at  moorings  M 1  and  M2. 

[26]  The  PAR  observations  at  moorings  (Figures  3a  and 
4a)  show'  that  the  incoming  solar  radiation  was  the  strongest 
during  the  extended  up  welling  event  (15-21  August).  This 


local  heating  interacted  with  cold,  upwelted  water  from  the 
nonhem  upwelling  center  at  Ano  Nuevo  [Ramp  et  al. ,  2005, 
Figure  1],  Little  cooling  is  observed  at  M2  (Figure  4b), 
which  is  located  in  an  area  only  weakly  affected  by  the  cold, 
upwelled  water  from  the  north  (mostly  moving  along  the 
entrance  to  and  into  the  Bay  between  M2  and  Ml).  Even  at 
Ml,  the  cooling  at  the  surface  takes  place  after  18  August 
when  there  is  a  reduction  in  the  PAR  (Figure  3a).  On 
average,  PAR  during  the  first,  brief  relaxation  event  (22— 
23  August)  is  lower  than  during  the  upwrelling,  especially  at 
the  mooring  M2  (Figure  4a).  However,  the  strong  increase  of 
the  daily  SWR  at  the  end  of  the  upwelling  and  the  beginning  of 
the  relaxation  events  (Figures  3a  and  4a)  leads  to  the  increase 
of  the  surface  temperatures  at  M 1  and  M2  during  this  brief 
period  of  relaxation  (Figures  3b  and  4b),  This  is  followed  by 
cooling  during  the  brief  upwelling  event  of  24-27  August. 
All  three  NCOM  ICON  runs  were  able  to  reproduce  the  above- 
described  sequence  of  cooling  and  warming,  especially  at 
mooring  M2  (Figure  4b).  The  model  simulations  during  the 
extended  relaxation  period  (28  August -2  September)  will  be 
discussed  in  detail  in  section  5.3. 

[27]  Statistics  for  runs  1-3  and  for  the  observations  at 
moorings  MI  and  M2  are  presented  in  Table  2.  Note  that 
there  are  potential  problems  with  the  use  of  point-to-point 
comparisons  with  mooring  observations  for  validation  of 
the  three  dimensional  circulation  models.  There  are  often 
spatial  and  temporal  shifts  between  model  predicted  and 
observed  features.  As  a  result,  relatively  small  spatial  offsets 
in  the  modeled  features  can  lead  to  an  unrealistically  poor 
assessment  of  model  predictions  when  comparisons  are 
done  with  the  point  mooring  observations, 

[2*]  All  three  runs  show  a  high,  larger  than  0.7  correlation 
with  observations.  Following  the  discussion  in  section  4,  the 
correlation  0.7  is  significant  if  the  number  of  degrees  of 
freedom  is  larger  than  or  equal  to  5,  AH  three  runs  have 
minimum,  maximum,  and  mean  values  of  temperature 
wanner  than  that  observed  at  Ml  and  M2.  This  is  a  result 
of  the  overestimation  of  SWR  fluxes  in  COAMPS  predic- 


Table  2,  Surface  Temperature  Statistics  15  August- 2  September 


Ml 

M2 

Min 

Max 

Mean 

SD 

Correlation 

Min 

Max 

Mean 

SD 

Correlation 

OBS 

Ii.ll 

15.02 

12.54 

0.84 

11,46 

16.62 

13,46 

1.36 

Runl 

12.60 

19.34 

14.69 

L73 

0.88 

14.10 

19.14 

14.77 

1.34 

0.81 

Run2 

14.22 

19,62 

16.00 

1.32 

0.87 

12.65 

20.31 

15.85 

1.81 

079 

Run3 

13.87 

16.27 

14,71 

0.50 

0.71 

13.16 

17.18 

15,03 

L.ll 

0,79 

Run4 

12.19 

17,54 

13.79 

1.29 

0.90 

13.91 

17.82 

15.30 

094 

0,73 

Run5 

12,16 

17J6 

13,68 

1.22 

0.89 

13.77 

17.40 

15.30 

0.87 

0,69 

Runf> 

13.60 

16,57 

14,65 

0.81 

0.82 

12,88 

17.58 

14.91 

1.06 

0,69 

Run? 

12,12 

16,53 

13,36 

1  IS 

0.82 

13,84 

16.74 

14.98 

0.82 

0,63 

Run8 

13.17 

16.09 

14,25 

0.66 

0.64 

13,33 

17.61 

15.18 

0,87 

0.73 
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Figure  5,  Observed  and  NCOM-1CON  model  predicted  mixed  layer  depth  (MLD)  comparisons  at  Ml 
and  M2. 


tions.  For  example,  in  run  I,  the  mean  temperatures  are 
higher  than  the  observed  means  by  2°C  at  Ml,  and  by  about 
l°C  at  M2.  The  maximum  and  minimum  values  deviate 
even  more  strongly  from  the  observed  values.  This  directly 
correlates  with  the  positive  biases  in  the  PAR  values 
presented  in  Figures  3a  and  4a,  The  mean  PAR  value  for 
the  “3km,  Cv3"  is  around  118  W/m2  at  Ml,  which  is  about 
14  times  larger  than  the  observed  mean  value  of  84  W/m2. 
In  order  to  correct  the  temperature  bias,  two  new  runs  were 
conducted  in  wrhich  the  values  of  the  SWR  flux  were 
adjusted  to  better  match  the  observed  values. 

[29]  Run  4  {Table  1 )  is  a  clone  of  run  1 ,  but  the  COAMPS 
SWR  fields  were  reduced  by  75  W/m2,  which  corresponds 
to  a  reduction  in  PAR  of  34  W/m2,  which  is  equal  to  the 
difference  between  the  CO  AM  PS-derived  and  observed 
values  of  PAR  at  Ml.  Run  5  (Table  t )  is  a  done  of  run  1, 
but  the  COAMPS  SWR  fields  w'ere  reduced  by  a  factor  of 
1.4.  In  both  runs,  the  mean  of  the  SWR  at  Ml  is  near  the 
observed  value  at  M  L  However,  reduction  by  a  factor  of  I A 
will  reduce  the  mean  and  variance  in  the  SWR  fields  in 
run  5,  while  in  run  4,  only  the  mean  value  is  reduced.  Note 
that  mean  values  of  SWR  in  runs  4  and  5  are  still  higher 
than  the  observed  values  at  M2,  The  statistics  from  runs  4 
and  5  are  presented  in  Table  2.  The  reduction  in  the 
COAM PS-predicted  SWR  significantly  reduced  the  mean 
and  maximum  values  of  the  NCOM  model  temperatures, 
while  maintaining  high  correlation  w  ith  the  observed  tem¬ 
peratures  at  Ml  and  M2.  Based  on  the  mooring  observa¬ 
tions,  correction  of  the  SWR  forcing  improves  agreement 
between  the  observed  and  the  NCOM  ICON  model-pre¬ 
dicted  SSTs. 

5*2.  Sensitivity  to  the  Vertical  Coordinate  System 

[30]  In  this  section,  we  compare  run  5  and  three  more  runs 
(Table  1 ).  Run  6  is  a  clone  of  run  2  but  is  forced  with  adjusted 
(divided  by  14)  SWR  from  the  “9km,Cv3*\  Run  7  is  a 
clone  of  run  5  (forced  with  “3km, Cv3”  forcing),  but  the 
NCOM  ICON  model  uses  a  hybrid  coordinate  system  in  the 
vertical  {see  section  2,  the  hybrid  grid  has  40  total  vertical 
layers  with  19  sigma  layers  and  2]  z-levels).  Run  8  is  a  clone 


of  run  7  but  is  forced  with  it9km,Cv3”  forcing.  In  both  runs  7 
and  8,  SWR  was  adjusted  as  in  runs  5  and  6+  Therefore, 
comparisons  of  runs  5-8  will  provide  insight  into  the 
sensitivity  of  the  model  results  to  the  vertical  coordinate 
systems  (sigma  versus  hybrid)  under  different  resolution 
atmospheric  forcing  (“3km,Cv3”  versus  “9km,Cv3"), 

[31]  SSTs  statistics  for  runs  5-8  (presented  in  Table  2) 
are  very  similar  with  slightly  better  correlations  for  runs  5 
and  6.  Figure  5  provides  comparisons  of  mixed  layer  depth 
(MLD)  for  runs  5-8,  MLD  is  defined  here  as  the  depth  at 
which  the  water  temperature  becomes  0.8°  less  than  the  SST 
[Kara  et  at,  2000].  At  MI,  run  5  reproduced  better  the 
observed  deepen ing/sh allowing  of  the  MLD  during  the 
upwell  mg/relaxation  events  than  the  other  considered  runs. 
This  is  reflected  in  the  statistics  presented  in  Table  3,  where 
the  correlations  with  the  observed  MLD  and  the  mot  mean 
square  (RMS)  errors  are  showm  for  runs  5-8.  The  MLD  for 
run  7  (“3km1Cv3”  forcing,  hybrid  coordinate  system)  shows 
a  noisy,  oscillating  MLD  during  strong  upwell  mg  (15-21 
August,  Figure  5).  Run  6  (“9km, Cv3“  forcing,  sigma  coor¬ 
dinates)  missed  the  deepening  of  the  MLD  during  this  event 
(Figure  5).  At  M2,  all  the  runs  show  a  deeper  than  observed 
MLD  during  upwelling  (15-21  August,  Figure  5),  This 
strong  deepening  in  the  model  results  might  be  an  indication 
that,  in  the  model  predictions,  upwelled  water  moves  off¬ 
shore  too  tar  by  reaching  the  location  of  M2.  Mixed  results  in 
the  prediction  of  M  LD  at  M2  are  reflected  in  Table  3,  where 
none  of  the  runs  shows  bener  results, 

5,3,  Simulations  During  the  Extended  Relaxation 
Event 

[32]  Both  runs  l  and  2  (forced  with  the  new'  COAM  PS 
“3km, Cv3”  and  “9km, Cv3"  products)  were  able  to  repro¬ 
duce  the  extreme  warming  during  the  extended  relaxation 
event  of  28  August-2  September,  while  run  3  (forced  with 
the  “9km,Cv2“  run)  missed  this  event  (Figures  3b  and  4b), 
One  of  the  reasons  for  this  is  that  the  atmospheric  forcing 
“9km,Cv2“  does  not  have  the  observed  increase  in  SWR 
during  the  relaxation  event  (Figures  3a  and  4a).  The 
“9km,Cv2”  was  also  used  to  force  the  POM  ICON  run 
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Table  3.  Comparisons  of  MLD  Predictions 


RMS  (m) 

Correlation 

Ml 

M2 

Ml 

M2 

RUN  5 

5.5042 

10,0535 

0,79 

0.54 

RUN  6 

6.0090 

1 1 ,885 1 

0,52 

0.64 

RUN  7 

10.7499 

12.8919 

0,54 

0.49 

RUN  8 

6.4523 

14.7254 

0,55 

0,56 

described  by  Ramp  et  al.  [2005]  and,  which  also  missed  the 
strong  temperature  increase  during  the  relaxation, 

[n]  As  was  speculated  by  Ramp  et  at ,  [2005],  the 
offshore  thermal  maximum  during  the  relaxation  period  at 
M2  is  ascribed  to  a  combination  of  advection  and  surface 
heating,  but  primarily  advection,  while  within  the  inner  bay 
(at  Ml),  the  thermal  maximum  was  due  to  surface  heating. 
This  corresponds  to  Rosenfeld  et  al.  [1994,  page  954], 
where  it  is  noted  that  a  large  surface  temperature  increase 
during  the  relaxation  “cannot  be  accounted  for  by  surface 
heat  flux  alone  and  must  be  result  of  advection”. 

[34]  Figure  6  shows  the  two-dimensional  structure  of  the 
wind  stresses  for  the  three  COAMPS  products  averaged 
over  the  relaxation  period.  It  is  dear  that  the  “9km,  Cv2” 
run  has  a  different  structure  than  the  “3km,Cv3”  and 
“9kmfCv3”  runs.  As  a  result,  the  offshore  California 
current  warm  water  was  advected  more  onshore  during 
the  relaxation  in  runs  l  and  2  with  the  “3km,Cv3”  and 
“9km, Cv3“  forcing  w'hen  compared  to  run  3  with  the 
“9km,Cv2“  forcing  (Figure  7).  This  is  another  reason 
why  runs  l  and  2  reproduce  the  strong  temperature  increase 
and  run  3  does  not.  Note  that  the  “3km, Cv3”  and 
“9km,Cv31'  wind  velocities  agree  better  with  the  aircraft- 
observed  wind  velocity  [see  Figure  8  in  Ramp  et  al ,  2005], 

[35]  Heat  contents  changes  for  runs  1  -3  were  estimated 
during  the  observed  surface  heating  event  of  30-31  August. 
The  area  of  the  bay  bounded  on  the  west  by  the  section 
crossing  the  location  of  the  mooring  Ml  (Figure  l),  and 
going  along  the  entrance  to  the  bay  from  Santa  Cruz  at  the 
north  to  Monterey  in  the  south,  was  considered  in  the  heat 


analysis.  During  30-31  August,  changes  in  the  heat  content 
due  to  SWR  for  NCOM  runs  1  and  2  (forced  with  the 
COAMPSv3)  are  about  30%  larger  than  run  3  (forced  with 
CQAMPSv2),  Note,  that  flux  due  to  SWR  is  around 
3.5  times  larger  than  the  total  flux  of  the  other  three:  latent, 
sensible,  and  long  wave  radiation  heat  fluxes.  During  30- 
3 1  August,  the  change  of  heat  content  due  to  advection  (for 
the  top  15  m)  is  about  3.5  times  larger  for  run  1,  and  about 
2  times  larger  for  run  2  in  comparison  to  run  3,  This 
supports  above  findings,  in  that  the  missed  increase  in 
SWR  and  the  deficiency  in  wind  stress  structure  for  the 
COAMPSv2  product  resulted  in  the  inability  of  the  NCOM 
ICON  run  3  and  the  POM  ICON  run  to  reproduce  the 
observed  extreme  local  heating  during  the  relaxation  period, 

5,4.  Simulations  During  the  Extended  Upwelling  Event 

[35]  Figure  8  presents  wind  stresses  averaged  over  the 
upwelling  period  (15-21  August)  for  the  “3km,Cv3“  and 
“9km,Cv3”  products.  The  corresponding  averaged  circu¬ 
lation  patterns  for  run  5  (forced  with  “3km,Cv3“)  and  run  6 
(forced  with  “9km,Cv3")  are  shown  in  Figure  9.  Both  runs 
reproduced  the  observed  upwelling  at  Ft.  Aho  Nuevo  and 
Pt.  Sur  (Figure  9),  At  the  same  time,  the  “3km,Cv3”  wind 
stress  is  stronger  than  the  9-km  stress  near  the  coast 
(Figure  8),  This  generates  stronger  upwelling  in  run  5  than 
in  run  6.  The  southward  jet,  between  the  wanner  water  of 
the  California  Current  and  the  upwelled  water,  flows  along 
the  entrance  to  the  bay  in  run  5.  This  agrees  with 
observations  from  Ramp  et  al  [2005]  and  the  CODAR- 
observed  surface  currents  presented  in  Figure  10  (left 
panel).  In  run  6,  the  surface  current  flows  offshore  and 
then  flows  along  the  bay. 

[37]  Run  6  shows  a  northward  flow  along  the  coast  that 
emanates  from  the  southern  open  boundary  of  the  model 
domain.  This  northward  flow  is  not  present  in  run  5.  Note 
that  during  the  relaxation  period  (Figure  7),  a  northward 
flow  at  the  surface  is  present  in  both  runs  1  and  2  (and  in 
runs  5  and  6,  not  shown  here).  During  relaxation,  when  the 
wind  forcing  is  weakening,  the  northward  flow  develops 
and  this  flow  is  well  observed  [Ramp  et  aly  2005],  The 


MEAN  WIND  STRESS  DURING  RELAXATION 
28  Aug  -  2  Sep,  2000 


123W  122W  123W  122W  123W  122W 


3km,  Cv3  9km,  Cv3  9km,  Cv2 

Figure  6.  Averaged  over  the  relaxation  event  wind  stresses  from  three  COAMPS  products:  (a) 
— 3km, Cv3;  (b)  -9km,Cv3,  and  (c)  — 9km,Cv2.  See  captions  of  Figure  2  for  explanations  of  notations. 
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Figure  7.  The  NCOM  ICON  model  predicted  SSTs  (averaged  over  the  relaxation  event)  from  runs  1  -3 
(Table  1)  forced  with  corresponding  COAMPS  wind  stresses  from  the  Figure  6, 


question  is  whether  the  surface  northward  flow  shown  in 
run  6  during  the  upwdling  (strong  winds)  is  a  physical 
phenomenon,  especially  since  run  5,  which  differs  only  in 
having  higher-resolution  atmospheric  forcing,  does  not 
reproduce  this  aspect. 

[js]  The  CODAR  radial  currents  observed  at  Pt+  Sur  and 
averaged  over  the  same  upwelling  event  (Figure  10,  right 
panel)  show  a  southward  flow  at  Pt  Sur,  which  corresponds 
to  the  results  of  run  5  but  not  run  6.  Therefore,  we  should 
conclude  that  the  development  of  northward  flow  at  the 
surface  during  the  upwelling  in  run  6  is  unrealistic.  What  is 
the  reason  for  the  development  of  this  artificial  northward 
flow?  Figure  1 1  shows  the  NCOM  CCS  velocity  normal  to 
the  southern  open  boundary  of  the  NCOM  ICON  domain. 


The  NCOM  CCS  velocity  is  averaged  over  the  upwdling 
event  Positive  values  indicate  northward  flow,  and  negative 
values  indicate  southward  flow.  Figure  1 1  demonstrates  that 
the  NCOM  CCS  velocity  has  southward  flow  at  the  surface 
close  to  the  coast.  Therefore,  open  boundary  conditions  at 
the  southern  boundary  of  the  NCOM  ICON  model  are  not 
responsible  for  the  development  of  the  northward  flow  in 
the  NCOM  ICON  run  6. 

[39]  Comparison  of  the  wind  patterns  of  the  “3km,Cv3" 
and  “9km, Cv3”  products  (Figure  8)  shows  the  weakening 
of  the  "9km,Cv3”  wind  stress  along  the  southern  coast  of 
the  NCOM  ICON  model  domain.  As  indicated  by  Capet  et 
al  [2004],  this  weakening  or  drop  off  at  the  coast  is  a  source 
for  intensification  of  the  northward  flow  south  to  Pt.  Sur  (as 


123W  122W  123W  I22W 

3km,  Cv3  9km, Cv3 

Figure  8.  Averaged  over  the  upwdling  event  wind  stresses  from  two  COAMPS  products:  (a)  -3km, Cv3 
and  (b)  -9km,Cv3. 
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MEAN  SST  AND  CURRENT 
UP  WELLING  15  Aug  -  21  Aug,  2000 
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Figure  9,  The  NCOM  ICON  model  predicted  SSTs  (averaged  over  the  upwelling  event)  from  runs  5 
and  6  (Table  1)  forced  with  corresponding  COAMPS  wind  stresses  from  the  Figure  8. 


it  is  seen  in  run  6).  Therefore  the  reason  for  the  development 
of  the  artificial  northward  flow  in  run  6  is  the  weakening  of 
the  9km  CO  AMPS  wind  at  the  coast  in  the  southern  portion 
of  the  domain.  The  impact  of  high-resolution  wind  fields  on 
the  POM  ICON  model  results  were  studied  by  Blencoe 
[2001], 

6.  Discussions  and  Conclusions 

[40]  Results  from  the  NCOM-based  model  (NCOM 
ICON)  of  the  Monterey  Bay  area  are  evaluated  during  15 
August-2  September  2000.  The  time  frame  includes  an 
extended  upwelling  event  during  15-21  August  and  an 


extended  wind  relaxation  event  from  28  August  through 
2  September. 

[41]  The  atmospheric  fields  from  two  COAMPS  versions 
(the  older  version  COAMPSv2  and  the  newer  version 
CO  AM  PS  v3)  are  compared  to  observed  wind  and  PAR 
values  at  moorings  MI  and  M2.  COAMPS  predicted  winds 
show  a  good  agreement  with  observations  in  the  sequence 
and  extent  of  each  up w el  ling -relaxation  event  observed 
during  15  August -2  September  2000.  As  a  result,  the 
NCOM  ICON  runs  forced  with  both  COAMPS  versions 
were  able  to  reproduce  the  sequence  of  the  observed 
upwell  mg/relaxation  events. 

[42]  However,  comparison  of  the  observed  and  COAM  PS 
estimated  PAR  values  show  that  there  is  overestimation  of 


UPWELLING  15  Aug  -  21  Aug,  2000 

CODAR  PREDICTED  CODAR  RAD1ALS 


Figure  10.  Averaged  over  the  upwelling  event  CODAR  derived  surface  currents  (left)  and  radials  at  Pt. 
Sur  (right). 
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MEAN  VELOCITY  AT  SOUTHERN  BOUNDARY 
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Figure  II.  The  NCOM  CCS  alongshore  velocity  plotted  along  the  southern  open  boundary  of  the 
NCOM  ICON  model.  The  velocity  is  averaged  over  the  up  welling  event.  Positive  values  mean  northward 
flow* 


SWR  (by  around  40%)  in  the  COAMPS  predictions.  The 
excessive  SWR  is  likely  related  to  a  modeling  of  low-level 
clouds,  which  are  very  extensive  during  summer  time  in  the 
Monterey  Bay  area*  As  a  result  of  the  overestimation  of 
SWR,  the  NCOM  ICON  model  results  show  minimum, 
maximum,  and  mean  values  of  temperature  warmer  than 
that  observed  at  Ml  and  M2.  In  order  to  conect  the 
temperature  bias,  the  NCOM  ICON  runs  were  conducted 
with  the  values  of  the  COAMPS  SWR  flux  adjusted  to 
better  match  the  observed  values.  It  is  shown  that  the 
correction  in  the  COAMPS-predicted  SW'R  significantly 
reduced  the  mean  and  maximum  values  of  the  NCOM 
model  temperatures,  while  maintaining  high  correlation 
with  the  observations  at  Ml  and  M2.  These  results  demon¬ 
strate  that  the  correction  of  the  atmospheric  model  SWR 
based  on  the  mooring  observations  improves  the  agreement 
between  the  observed  and  the  NCOM  predicted  SSTs  at  the 
mooring  locations. 

[43]  The  NCOM  ICON  runs  forced  with  different  reso¬ 
lution  atmospheric  forcing  (3-  versus  9-km  resolution)  as 
well  as  the  NCOM  ICON  runs  employing  different  vertical 
coordinate  systems  (sigma  versus  hybrid)  do  not  show 
significant  differences  in  the  prediction  of  SST  (Table  2)* 
However,  comparisons  of  mixed  layer  predictions  show> 
that  the  run  5,  with  a  sigma  vertical  coordinate  system 
and  forced  with  3 -km  resolution  forcing,  better  reproduced 
the  observed  deepening/shal lowing  of  the  MLD  during 
the  upwelling/relaxation  at  the  MI  mooring  location.  At 
mooring  M2,  all  the  runs  show  a  deeper  than  observed 
MLD  during  upwelling  (15-21  August).  This  strong  deep¬ 
ening  might  be  an  indication  that,  in  the  model  results, 
upwelled  water  moves  offshore  too  far  by  reaching  the 
location  of  M2.  Mixed  results  in  the  prediction  of  MLD  at 
M2  are  reflected  in  Table  3,  where  none  of  the  runs  show 
better  results. 

[44]  During  the  wind  relaxation  event  (28  August- 
2  September),  strong  onshore  translation  of  warmer 
California  Current  water  combined  with  local  surface  heat¬ 
ing  inside  the  bay  resulted  in  extraordinary  surface  wanning 
within  the  bay.  As  reported  by  Ramp  et  ai  [2005],  the  POM 
ICON  model  simulations  forced  with  the  older  version  of 
the  atmospheric  model  COAMPS  (COAMPSv2)  were  not 


able  to  reproduce  this  observed  extraordinary  surface  warm¬ 
ing.  Similarly,  the  NCOM  ICON  run  forced  with  the  same 
version  of  the  COAMPS  missed  the  event.  Analysis  of  the 
CGAMPSv2  atmospheric  fields  show  two  reasons  for  this: 
the  missed  observed  increase  in  SWR  during  the  relaxation 
event  and  the  deficiency  in  wind  stress  structure,  which 
does  not  bring  offshore  warmer  water  into  the  Bay.  How¬ 
ever,  the  NCOM  ICON  model  runs  forced  with  the  new 
COAMPSv3  fields  (3-  and  9-km  resolution)  were  able  to 
reproduce  the  observed  event. 

[45]  During  the  extended  upwelling  event,  the  NCOM 
ICON  runs  forced  with  3-  and  9-km  resolution  COAMPSv3 
fields  (ams  5  and  6,  Table  l)t  respectively,  showr  differences 
in  surface  circulation  patterns.  The  most  distinct  difference 
is  in  the  southern  portion  of  the  NCOM  ICON  model 
domain:  the  run  forced  with  9-km  forcing  shows  the 
northward  flow  along  the  coast  that  emanates  from  the 
southern  open  boundary  of  the  model  domain.  This  north¬ 
ward  flow  is  not  present  in  the  run  forced  with  3-km 
forcing.  Observations  like  CODAR  radials,  as  well  as 
predictions  from  the  larger  scale  California  Current  model 
(NCOM  CCS  predictions  which  are  used  on  open  bound¬ 
aries  of  the  NCOM  ICON  model),  show  southward  flow  in 
the  southern  portion  of  the  model  domain.  Comparison  of 
the  wind  patterns  of  the  3-  and  9-km  products  (Figure  8) 
shows  the  weakening  of  the  9-km  wind  stress  along  the 
southern  coast  of  the  NCOM  ICON  model.  This  weakening 
is  responsible  for  the  development  of  the  artificial  north¬ 
ward  flow  in  the  NCOM  ICON  run.  The  weakening  of  the 
9-km  wind  stress  at  the  coast  is  also  indicated  as  a  source  tor 
intensification  of  northward  flow  in  the  work  of  Capet  et  al 
[2004]. 
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